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Abstract 

We present a refined and improved study of the influence of screening on the effective fine 
structure constant of graphene, a *, as measured in graphite using inelastic x-ray scattering. This 
follow-up to our previous study 1 ! was carried out with two times better energy resolution, five 
times better momentum resolution, and improved experimental setup with lower background. We 
compare our results to RPA calculations and evaluate the relative importance of interlayer hopping, 
excitonic corrections, and screening from high energy excitations involving the a bands. We find 
that the static, limiting value of a* falls in the range 0.25 to 0.35, which is higher than our previous 
result of 0. l4D, but still below the value expected from RPA. We show the reduced value is not a 
consequence of interlayer hopping effects, which were ignored in our previous analyst, but of a 
combination of excitonic effects in the n —> ir* particle-hole continuum, and background screening 
from the u-bonded electrons. We find that cr-band screening is extremely strong at distances of 
less than a few nm, and should be highly effective at screening out short-distance, Hubbard-like 
interactions in graphene as well as other carbon allotropes. 
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I. INTRODUCTION 


Since its synthesis a decade ago, graphene has been the subject of intense research ! 2 3 . 
In the long-wavelength limit, graphene satisfies the massless Dirac equation, where its elec¬ 
trons follow a linear band dispersion given by E = ±hvfq , where the ± refers to the con¬ 
duction/valence bands, Vf ~ 10 8 cm/s ~ c/300 is the Fermi velocity, and q is the electron 
momentum (in 2 D)HGD As a result of this dispersion, the bare electron-electron interaction 
strength - as defined by the fine-structure constant, a = U/K & 2.2/e, where e is the 
dielectric constant of the substrate on which the graphene is mounted (a ~ 2.2 for sus¬ 
pended graphene) - is of order unity, in contrast to the interaction strength in QED, where 
a = 1/137; thus, many-body effects should, in principle, play a significant role in graphene 6 . 
In particular, the screened value of a, which one could call a renormalized or effective fine 
structure constant, a*, is a subject of great theoretical interestP®, and is relevant to all 
Dirac systems including topological insulator surface states®^, some classes of transition 
metal dichalcogenide^, three-dimensional Dirac and Weyl semimetals- 415 , etc. 


Experiment and theory, however, give conflicting views on the value of a*. Using inelastic 
X-ray scattering (IXS) experiments performed on graphit^, we previously found that, for 
freestanding graphene, a* ~ 0.15 as q —> 0 and oo —> 0, and cited excitonic shifts in the 
7 r —y tv* , particle-hole continuum as the origin of the reduced value of a*. Interpretation 
of this experiment was based on the assumption that the interaction between graphene 
sheets in graphite is primarily Coulombic, with interlayer hopping playing a secondary role 1 . 
The experiment, however, had an energy resolution of only 0.3 eV, when the features of 
primary interest range from 1-15 eV. Additionally, the lowest momenta measured were 
at 0.238 A ^ 1 and 0.476 A -1 so the extrapolation of q —> 0 relied on only a few points. 
Subsequent calculations in the random-phase-approximation 16 (RPA) using the full 7 r-bands 


and including the interlayer hopping parameters y 3 and 71 in graphite (see Fig. 2a) suggested 
that this discrepancy was due to graphitic effects, i.e., the presence of interlayer hopping 
effects ignored in our analysis. 


In this paper, we present a refined measurement of a* in graphene as measured in graphite, 
using an improved version of the previous methocP, with higher energy and momentum 
resolution, as well as reduced elastic background. We evaluate the importance of interlayer 
hopping, excitonic corrections, and the higher energy cr-bands to the effective fine-structure 
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constant. We compare to theoretical results obtained within the RPA, and show that for 
q ~ 0.212 A -1 (length scale ~ 30 A), a* falls in the range 0.25 to 0.35, which is well 
below the RPA value for graphene. We show the origin of this anomalously low value is not 
graphitic (i.e., interlayer hopping) effects, but a combination of screening from the cr-bonded 
electrons, which contribute significantly at finite momentum, and excitonic effects in the n 
bands, which comprise a beyond-PRA correction to the screening. While the effect of a 
band screening decreases with decreasing q, we argue such effects remain relevant to length 
scales up to rs-/ 30A. 


II. EXPERIMENT 


IXS measurements were carried out in Sector 9 at the Advanced Photon Source. Energy 
analysis was done with a diced, Si-444, spherical backscattering analyzer operating at 7.81 
keV, which imaged scattered photons onto a MYTHEN microstrip detector. The total 
energy resolution of the instrument was 175 meV, which is a factor of two better than our 
previous studj^. 

For these measurements, we constructed a new sample chamber equipped with a moving, 
in-vacuum beam stop, providing access to scattering angles as low as 1°. The momentum 
resolution was set to 0.15-0.3 A -1 , depending upon the scan region. The chamber was 
designed for windowless operation with the synchrotron beam pipes, reducing scattering 
from upstream vacuum windows. The final setup exhibited greatly reduced background in 
the 1-5 eV range, particularly from the quasielastic line, a key source of error as described 
in Appendix [A] 

Experiments were done on both ZYA-grade, highly-ordered pyrolitic graphite (HOPG) 
and high quality, commercially obtained graphite single-crystals. The spectra of the two 
were found to be indistinguishable in the momentum range in this study. The experimental 
data were placed in absolute units of (eV A 3 )^ 1 by normalizing using the F sum rule, as 
described in Ref. 1 and Appendix [A] 
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III. CONVERSION OF EXPERIMENTAL DATA 


In order to obtain information about a *, we need to determine the two-dimensional, 
density-density response function of graphene, X 2 d{q,w), which is the most general descrip¬ 
tion of the collective charge dynamics of the systenP. IXS measurements cannot be carried 
out on single-layer graphene, however, so a means of determining properties of graphene 
from measurements on graphite is required. 

The measured intensity in an IXS experiment at a given momentum and energy is pro¬ 
portional to the dynamic structure factor, S(q,cu), which is related to the density-density 
response function of three-dimensional graphite, X 3 d{q,w), by the quantum mechanical ver¬ 
sion of the fluctuation-dissipation theorem, 


S(q,u) = -- 


■Imx3£>(g,w). 


7r 1 — e~ 

We note that X 3 D is causal and satisfies the Kramers-Kronig relation 

Imx 3 D(q,oS) 


2 f c 

Rex 3 D (q,u) = -V \ 

71 Jo 


du 


( 1 ) 


( 2 ) 


(c o ') 2 — U) 2 

So, by extrapolating to u> —> oo, interpolating the discrete points in to, and integrating, it is 
possible to determine the full, complex X 3 D(q,w) for graphite in absolute units of (eVA 3 ) -1 , 


as shown in Fig. la Note that the onset of tt —>• n* particle-hole excitations is visible in 
the spectra near uj = Vfq, indicating that we are probing the relevant, low-energy valence 
excitations. The peak at u ~ 8 eV is the well-known 7r-plasmon, which is not a free carrier 
plasmon but an effect of the Van Hove singularity at the top of the n band 1 . 

Learning about graphene requires hireling a relationship between its density response, 
X 2 n(q, to), and the measured, three-dimensional response of graphite, X 3 D(q, oj). This can 
be accomplished by making three assumptions: (1) that the primary interaction between 
layers in graphite is electrostatic, with interlayer hopping playing a negligible role, (2) that 
the layers are arbitrarily thin, and (3) that the electrons are distributed homogeneously 
within each layer. In this case, the two response functions are related by 1 


/ N_ X3 D(q,u)d 

X2D 1 _ - F(q)}x3D(q,u)d' 

Here, V 2 o(q) = 27re 2 /g is the two-dimensional Coulomb interaction, and F(q) = sinh(qd)/(cosh(qd) — 
cos (q z d)) is a structure factor that describes the Coulomb interaction in a layered, three- 
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dimensional system, i.e., V 3 d{q) = V 2 D(o)F(q)d, where d = 3.35A is the interlayer 
spacing 1 ® 7 i 

Eq. [3] makes use of the fact that, if the above three assumptions are valid, the polarization 
function of graphene, II 2 d(< 1 ,w), is (apart from units) identical to that for graphite,^ i.e., 


n 2 D (q,u) = U 3D (q,uj)d. 


( 4 ) 


These quantities are related to the respective susceptibilities by Xsd = n 3 d/^ 3 d and X‘ 2 D = 
^■ 2 d/^ 2 D, where e 3 o = 1 — V 3 oIl 3 o and 62 d = 1 — are the dielectric functions. 


The spectra for ImygD obtained in this manner are shown in Fig. [Tbj and the associated 
polarization functions in Fig. [3} Crucially, the X 2 d(q , w) obtained from Eq. [3] exhibits 
the correct asymptotic properties expected for two-dimensional graphene (see Fig. [4a] and 
Section VI). 

Written in terms of the other quantities, the screened, effective fine structure constant of 
graphene a* is 


(X 

a*{<hu) = -= a ■ [1 + V 2D (q)x 2 D(q,uj)] . (5) 

^2 D 

This quantity can be thought of as the value of the fine structure constant incorporating 
screening corrections to all orders in perturbation theory. Note that, although Eq. [5]appears 
similar to an RPA expression, we have at no point assumed that only RPA polarization 
bubbles contribute to the susceptibility; indeed, the x we recover from experiment includes 
all screening processes including excitonic effects and other corrections beyond the RPA. 


IV. RPA CALCULATIONS 

To aid interpretation of the experimental data, we performed a 7 r-electron tight binding 
calculation of the susceptibilities of graphite and graphene in the random phase approxima¬ 
tion (RPA), following the approach of Yuan et al. lh . In this approach, one calculates the 
polarization functions for single-layer graphene (SLG), IRd, and multilayer graphite (MLG), 
II 30 , using the Lindhard formula, as described in Appendix B. The polarization functions 
are then used to determine the dielectric functions and susceptibilities via the relationships 
discussed in Section III. 
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The RPA approach neglects excitonic effects, which have been argued to be crucial for 
interpreting the particle-hole excitation spectrum of graphene 118 . However, RPA has the 
advantage of allowing one to switch on and off the interlayer hopping parameters, 71 and 


73 defined in Fig. 2a, allowing one to evaluate the consequences of both interlayer hopping 
and excitonic effects by comparing to the experimental IXS data. 

Like Yuan 1 ®, our calculations are based on a tight binding theory that includes only the 
electrons in the 7 r bands (see Appendix B). The effects of the deeper a bands are incorporated 
phenomenologically as a frequency-independent, background dielectric constant 


Ac 0 + 1 - (/c 0 - l)e qLN 


( 6 ) 


k 0 + 1 + (k 0 - l)e _ljLj v 

where k 0 = 2.4 is the background dielectric constant of graphite 2 ^, L N = d m + (N — 1 )d is 
the height of a multilayer graphite system with N layers, and d m = 2.8 A is the thickness of 
a single graphene sheet 16 . Eq. [b] exhibits the proper scaling with the number of layers at 
low momenta, and was argued to correctly describe screening by excitations involving the 
sigma bandit 1 ®. Note that for infinite, MLG, this expression reduces to /too (q) = f° r any 
nonzero momentum. 

We differ from Yuaii 1 ®, however, in our use of Kjy(q). In Ref. 1 ^, Kjv(q) was used only to 
reduce the strength of the Coulomb interaction, by using a screened W(q) = V(q)/n(q). 
As shown in Appendix C, this usage neglects interference between screening by n and a 
electrons. The correct use of kN^ q) is in the relationship between the polarization and 
dielectric functions, 


e(q,u) = K N (q)-V(q)n(q,u), (7) 

where N is chosen to match the dimensionality of the system. 

Having determined the polarization functions for MLG, Usoiq,^), and SLG, Uo£)(q,uj), 
we calculated four different susceptibilities for comparison to experiment. The first, 

R pA ( \ _ n 3D (g,q;) _ 

3D ’ Yx>(<z) — V2D(q)F(q)dU 3D (q,uj) ’ 

is the three-dimensional susceptibility for MLG in the RPA, including the effects of both 
interlayer hopping and electrostatic coupling between the layers. The imaginary part of this 
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quantity could be compared directly to the experimental data to assess the importance of 
excitonic effects. The second, 


x2E A (q,<*>) = 


fl 2 d(<?,w) 


- (9) 

ki (q) - V 2 D(q)n 2 D(q,v) 

is the two-dimensional susceptibility for SLG. This quantity could be compared to the con¬ 
verted experimental data resulting from the application of Eq. 3 to assess the combined 
importance of excitonic and interlayer hopping effects. The third quantity 


xSS A (q,u) = 


U 3D (q,uj)d 


( 10 ) 


Ki{q) - V 2 D (q)II 3 D (q,cj)d 
is a two-dimensional propagator describing, in RPA, how the physics of graphene would 
be modified if 71.3 7 ^ 0 , i.e., if it were subjected to the same interlayer hopping effects as 
graphite. Finally, 


-hpa/ \ = _ II 2 D (q,u;)/d _ 

(9 ’ J Koc(q)-V 2 D (q)F(q)n 2 D (q,co)’ 

is a three-dimensional propagator for graphite in which the interlayer hopping is switched 
off, i.e., 71,3 = 0, but the layers are still coupled by electrostatic effects. By comparing 
these quantites to one another, to the experimental data, and to the converted data deter¬ 
mined from Eq. 3, it should be possible to disentangle the contribution of excitonic effects, 
interlayer hopping, and cr-band screening on the effective fine structure constant of graphene. 



V. RESULTS 

We are now in a position to evaluate the source of the anomalous screening observed 
in our original study 1 . In particular, we wish to know whether this screening is due to 
the graphitic effects discussed recently 1 ^, i.e. the interlayer hopping parameters 71 and 73 , 
background screening from the cr-bands, or excitonic effects as we claimed previously 1 . 

In Fig. we compare the experimental data to the imaginary part of \ 3 D computed 
from RPA both with and without interlayer hopping (Eqs. 8 and 11, respectively). The 
first thing to note is that RPA does a poor job at reproducing the energy of the so-called 
“ 7 r-plasmon” at ~ 8 eV. Note that this excitation is not a free carrier plasmon, but a result of 
the Van Hove singularity at the top of the n band 1 . We attempted to improve the agreement 
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by adjusting the parameters t, 71 , 3 , and /c 0 , but any parameters that produced the correct 
plasmon dispersion were dramatically different from the commonly accepted values 6 22 '' of 
t — 3 eV, 71 = 0.4 eV, 73 = 0.3 eV, and k 0 = 2.4. Changes to these parameters also 
drastically distorted the qualitative shape of the tt —» tt* continuum. The discrepancy is 
most likely due to excitonic effects not captured by RPA 1? . We therefore show only the 
spectra calculated using the commonly accepted values for the hopping parameters and k 0 . 
Note that, in agreement with Ref. 1 ^, the effect of a nonzero 7^3 is not to qualitatively change 
the spectra, but to smear some of the features near threshold. Such effects are subtle and 
not easily distinguished from broadening due to experimental resolution. 

In addition, a significant discrepancy is observed between the size of the experimental 
response and that computed with RPA (Fig. la). The latter is significantly larger in 
magnitude and had to be multiplied by an arbitrary scale factor, ranging from 0.15 to unity, 
to allow visual comparison with the experiment. This magnitude discrepancy is fundamental, 
and arises from the significant spectral weight in the cr-bands in the experimental data 
that is absent from the RPA result. Despite residing at ~ 25 eV, the cr-bands are quite 
polarizable and have significant influence on the value of x a t low energy, because x is a 
nonlinear function of II 22 23 . In other words, our RPA approach treats the system as if all 
the screening is accomplished by the n electrons, while in reality some of the work is done 
instead by the cr-bonded electrons. RPA therefore predicts that the response from the n 
electrons is anomalously large, since the added cr screening channel is absent. The influence 
of the cr-bands increases with increasing q, reaching a maximum at around 2 A~ l . Note 
that this problem is not corrected by including a background Kn(q), which does not replace 
spectral weight that is missing from the F sum rule (note that using inside a screened 

W(q), as was done in Ref. 1 -^, does not improve the situation). 

Still looking at Fig. la, another key observation is that the onset of the 7 r —y n* continuum 
in the experiment is consistently lower in energy than in the RPA calculation. These shifts 
are even better visible in the imaginary part of the polarization function (Fig. 3), which is 
better reflective of single particle excitations than the susceptibility, which emphasizes the 
collective modes 24 . For momenta q > 0.35A~ 4 , these shifts are as much as 1 eV compared 
to RPA, which is too large to be explained by graphitic effects, and must be attributed 
to excitonic interaction between the valence hole and photoexcited electron, as discussed 
previously 1 ? . At lower momenta, the shifts are smaller and similar in magnituide to what is 






expected from interlayer hopping effects, though excitonic effects surely still play some role. 

We turn now to the two-dimensional response, Imx 2 D{q,u), shown in Fig. lb, which 
compares the experimental data converted using Eq. 3 to the two-dimensional RPA response 
calculated with and without interlayer hopping (Eqs. 9 and 10). The discrepancy in the n 
plasmon energy observed in Fig. la, as well as the excitonic shifts in the n —> 7 r* continuum, 
are visible also in the two-dimensional response. However, the magnitude discrepancy that 
is so pronounced in three dimensions is much less significant in two dimensions and fades 
as q decreases. At the lowest momentum point, q = 0.212A -1 , the magnitude discrepancy 
is absent entirely. This means that in two dimensions, in contrast to the three dimensional 
case, the polarizability of the high energy o bands ceases to be important in the limit of 
low momentum. This confirms one’s intuition that screening from high energy excitations 
should be unimportant in two dimensions as q —> 0 . 


VI. DISCUSSION 

We now turn to the fundamental question of how the effects of interlayer hopping, exci¬ 
tonic shifts in the spectra, and the polarizability of the o bands influence the value of the 
screened, effective fine structure constant, a*(q,ou) = ao[l + V2D(q)X2D(q,u)\- The quantity 
a*(q, uj) is a strong function of both energy and momentum, and its value differs greatly 
depending upon the relative size of the energy, c o, and the momentum, qvF■ For transport 
experiments in a micron-sized device at low temperature, oj/vpq ~ ksTL/2'KfiVF ~ 10 -3 . 
So the quantity of primary interest is the zero-frequency quantity o.*(q, 0) in the limit of 
small momentum. We note that this quantity is purely real, and could be either larger or 
smaller than the bare value, a:o = e 2 /eo Hvf- 

In Fig. 4a, we plot the momentum dependence of the ratios D(q, 0 )/qd and X 3 D(q, 0 )/q 
for both the experimental data as well as the RPA calculations both with ( 7^3 ^ 0) and 
without ( 7 i )3 = 0) interlayer hopping. Note that both of these quantities are real and have 
units (eUA 2 ) -1 , allowing them to be compared on a single scatter plot. As expected, the 
quantity X 3 D(q, 0 )/q tends toward 0 as q —> 0 , since the susceptibility in three dimensions 
is expected to scale like q 2 at small q. All three X2D(q,0)/qd curves, on the other hand, 
converge to a constant as q —» 0 , indicating that the quantity X 2 D(q, 0 ) scales linearly in q 
at low momenta, which is expected in two dimensions. This scaling is crucially important 
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for two reasons. First, it demonstrates that our conversion expression, Eq. 3, gives the 
proper asymptotic behavior of the susceptibility in two dimensions, which is a powerful 
validation of our method. The second is that the product V 2 d(q) • X 2 d(q^) converges to 
a constant at low q for all momenta, indicating that graphene exhibits a finite screening 
strength, e(q, 0 ) = [1 + V2D(q)X2D(q, 0 )] -1 , in the long-wavelength limit. 

An important conclusion about interlayer hopping can immediately be reached from Fig. 
4a by comparing the RPA calculations of X 2 d{q, 0 )/qd with and without interlayer hopping. 
Notice that the curves are only distinguishable for momenta less than about 0.2 A -1 , which 
is below the lowest momentum measured in our experiment. What this implies is that, while 
graphitic effects can in principle influence the value of a* at low momenta, such effects are 
likely irrelevant in the regime of our experiment, and that excitonic corrections are in part 
responsible for deviations compared to RPA visible in Figs. la,b at momenta q < 0.35A -1 . 

Fig. 4a also reveals a major discrepancy between the experiment and theoretical curves at 
large momenta. For q > O.bA^ 1 , the RPA curves drop off gradually, while the experimental 
points increase dramatically. This increase is a consequence of screening from the a bands, 
which exhibits a broad maximum at q ~ 2 A _1 , which is of the order the carbon-carbon 
bond length 1 . This discrepancy is, again, a consequence of the absence of a band electrons 
in the RPA calculation. Hence, the quantity n(q) taken from Ref. 16 } while argued to be a 
good description at small q, fails to describe any of the major features of a band screening 
at larger momenta. Fortunately, the influence of these high energy excitations on \2D fades 
at small momentum, suggesting that we should still get useful information about the low -q 
value of a* from our study. 

Finally, in Fig. 4b we plot the magnitude of the static, effective fine structure constant, 
a*(q,0), as a function of momentum for the experiment (i.e., deduced from Eq. 3), and for 
RPA both with and without interlayer hopping. For reference, the points from our previous 
study 1 are also shown. The difference between the two experiments can be attributed to 
better resolution, better statistics, and lower experimental background in the current study. 
For momenta q > 0.5A -1 the value of a* is dominated by screening from the cr bands. 
Unfortunately, the error bars on the lowest momentum point, at which cr-band effects are 
negligible, are extremely large due to extrapolation of the F sum rule (see Appendix A). The 
remaining points are converging toward a value in the range of 0.25 to 0.35. This quantity is 
larger than given in our previous study, but still smaller than the value expected from RPA, 
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and implies an asymptotic dielectric constant in the range 6.2 to 8 . 8 . Because interlayer 
hopping effects are irrelevant in this regime, we conclude that the anomalously low value 
of a* is a combined consequence of a band screening and excitonic shifts in the 7 r —)• 7 r* 
continuum. 


VII. CONCLUSIONS 

We presented significantly improved measurements of the momentum- and energy- depen¬ 
dent effective fine structure constant of graphene, a*(q,u), using inelastic x-ray scattering 
measurements of graphite. We deduce an asymptotic value of a* in the range 0.25-0.35, 
which is larger than stated in our previous study 1 but smaller than the RPA value of 0.49, 
and implies an asymptotic value of e in the range 6 . 6 - 8 . 8 . We also performed 7 r-band RPA 
calculations with and without interlayer hopping, and found that graphitic effects from in¬ 
terlayer hopping have no significant effect on screening in the range of momenta we studied. 
Screening by a band electrons, on the other hand, contribute significantly at finite wave 
vector. We therefore conclude that the anomalously low value of a* may be attributed to a 
combination of excitonic effects in the 7 r —y n* continuum and a-band effects. The latter are 
very large at finite q, and should strongly screen short-ranged, Hubbard-likc interactions in 
graphene and carbon allotropes more generally. 

We close by commenting on the discrepancies between our conclusions and the RPA 
calculation described in Ref. 16 . That study reported a calculated static dielectric constant 
at q = 0.238A- 1 that was not far from our earlier result derived from IXS experiments on 
graphite 1 . They attributed the anomalously low value of a* to graphitic effects from the 
interlayer hopping. Here we have shown that the interlayer hopping parameters, 71 and 73 , 
have no effect on the dielectric response at the lowest momenta measured in our experiments. 
The discrepancy between these two conclusions arises from the way the Coulomb interaction 
was handled in the two studies, which differ in two important ways. First, Ref . 16 ignores 
the interaction between screening in the cr and 7r channels, which must be accounted for as 
described in our Appendix C. Second, in comparing to our result, they failed to eliminate 
the Coulomb interaction between the layers, which we have done through the use of Eq. 3. 
These differences result significantly larger values of e(q, 0) in Ref. 1 ^, which they atribute to 
interlayer hopping effects. What we have shown is that, if the Coulomb effects are handled 
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properly, the apparent effects of interlayer hopping are in fact small, and the primary causes 
of increased screening are excitonic effects and a band screening. 
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Appendix A: Data Processing 


Three steps of data processing are required to extract the response function, % 3 d(q,w), 
from the raw experimental spectra. The first is substraction of the quasiclastic line, which is 
very intense in the experiment. The second is scaling the data using the E-sum rule, which 
provides the imaginary part, lm.X 3 D(Q, cu), in absolute units. The third is Kramers-Kronig 
analysis, which is needed to determine the real part, Re X 3 d(q, w)- 

To accomplish the first of the three, we used a pseudovoigt function—a linear combination 
of a Gaussian and a Lorentzian—to model the quasielastic lineshape: 


( x ~ x o)^ . 

ce 2 a‘-> + (1 — c) 

\2 


y = b + Ax^ 

ce + (1 - c) ^__ , X >= X 0 . 


(x-xo ) 2 +l 2 ’ 
Vv) 2 


X < x 0 


(Al) 


Here, a is the Gaussian width, l = y / 21og2cr is the Lorentzian width, c is a relative amplitude 
of the Guassian and Lorentzian components, Xo is the center of the lineshape, and b and A are 
overall background and amplitude factors, respectively. The parameter r] is an asymmetry 
factor that allows the positive- and negative-energy sides of the lineshape to be different, 
which is a feature in the experiment due to energy-dependent background scattering from 
the sample chamber. In particular, the Lorentzian width l is defined this way to best 
approximate a voigt function, which is a convolution of a Gaussian and Lorentzian. In 
our experiments, background measurements on lithium fluoride, a large band gap material, 


showed b = 0 for our setup, so this quantity was held fixed in our fits. Figs, pal and 5c show 
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subtractions for the two different momentum transfers q = 0.212 A -1 and q = 0.282 A -1 , 
respectively. To estimate the error introduced by subtracting the elastic line, we repeated 
the procedure with different constraints and initial parameters and recorded the range of 
values achieved in optimized fits. 

To normalize our data, we applied the F-surri rule over the entire momentum range. The 


results are shown in Fig. la. Energy-loss spectra from the current study were measured only 
over the interval —4 < u < 15 eV, so it was necessary to supplement using data from Ref. 1 
to achieve a sufficiently wide range of energy to perform a Kramers-Kronig integration. This 
was complicated by the fact that the q values in the current study are not the same as those 
in from Ref. 1 . To overcome this mismatch, we linearly interpolated the data from Ref. ^ to 
match our q values. To properly normalize and patch our spectra to the previous spectra, 
we first note that in Ref. 1 ', the proportionality between the measured spectra and S(q, u) 
(and, thus, x(qr,u;)) was determined by fitting the spectra to the F-surri rule. In our new ex¬ 
periment, for each separate experimental run, we took our spectra, with the exception of the 
q = 0.212 A -1 , and fit to the spectra of Ref.^, finding a multiplicative factor for each momen¬ 
tum transfer. We then used the average of these multiplicative factors to scale the complete 
data set. This was done because the multiplicative factor is determined predominantly by 
incident beam features including intensity, angular divergence, and longitudinal bandwidth, 
as well as sample thickness, which are features that are constant or nearly constant within a 
single experimental run, but may vary from one beam time to the next. For q = 0.212 A -1 
we took the average multiplicative factor for the q = 0.352 A -1 and q = 0.563 A -1 spectra 
(the other two momentum transfers from that run) and defined that its multiplicative fac¬ 
tor rather than directly fitting to the old data; this is because the lowest momentum value 
from the old data was q = 0.238 A -1 , and extrapolation to lower momenta is not reliable. 


Figs. 5b and 5d show our final, normalized spectra for q = 0.212 A -1 and q = 0.282 A” 1 , 
respectively, in absolute units. After acquiring these spectra, Kramers-Kronig analysis was 
straight forward using standard methods. 

The main advantage of our data processing procedure is that it fits the 7r-plasmon inten¬ 
sities well, which is ideal given that 7r-plasmon is the strongest contribution to the Kramers- 
Kronig integral in the energy range we measure. Unfortunately, for a few momentum values, 
our method introduced discontinuities at the patch energy (« 15 eV), as illustrated in Fig. 


5cl There being no objective way to correct for this, we simply included the effect in our 
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estimate of the error bars on the value of x(q, w) (see Fig. la). Indeed, the spectra in Fig. [I] 
represent the average values of x after elastic line subtraction and intensity normalization, 
and the error bars are a sum in quadrature of the maximum positive or negative deviation 
from the mean and a 5% estimate of the variance of the sum-rule normalization. 


Appendix B: RPA 


As discussed in the main text, RPA calculations for graphite and graphene are based on 
a tight-binding model of the 7 r electrons. The graphite calculations assumed ABA(Bernal) 


stacking, illustrated in Fig. 2a including nearest-neighbor interlayer hopping. The tight- 


binding Hamiltonian has the explicit form 


N, 


layer 


H = ~ t J2 [ a h b B + b h a *>' 

l(*j> 


7i 


£ 

leven,{{j,j'}) 






b \,i a :',i+i + a ]',i+ibj,i 


73 

£ 

b],i a j',i+ 1 + a ]',z+i bj,i 


^even 



+ a ],fir,i+i + b ]’,i+i a j,i 


(Bl) 


where the sum is taken over nearest neighbors (i , j) and the layers, l, t — 3 eV is the 
in-plane hopping parameter, — 0.4 eV is the vertical hopping between the A and B 
sublattice (indicated by the double brackets), and 73 = 0.3 eV is the nearest A and B 
interlayer coupling (indicated by the triple brackets) (see Fig. 2a). Single layer graphene 


can be described by the first sum in Eq. Bl For graphene, the wavefunctions can be 
computed explicitly, and we can write the overlap function as 

4>*{k) (j>(k + q) 


f»A k ,q 11 ) = 


|0(fc)| \(j)(k + q] 


where 


+ s ■ s' 


'V3, 


(f>{k) = e~ ik * a + 2 e-^ a/2 cos ( -y k y a 


(B2) 


(B3) 
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and a = 1.42 A is the in-plane carbon-carbon distance. 

Once the tight-binding wave functions were obtained, the polarization functions were 
obtained by numerically evaluating the Linclhard formula, 


n(m|,cn) = 


9s 


(2tt) 


d 2 k 


IBZ 


X 


ST f (k- n \ n F[E s (k)} - n F [E s \k + q)} 
‘~_ 1 SS ’ " E s (k) — E s '(k + q) + u + irj' 


(B4) 


where the integral is over the Brillouin zone, d is the dimensionality (either 2 or 3 for graphene 
or graphite, respectively), g s is the spin degeneracy, E s / s '(k) is the energy dispersion for the 
conduction (s = 1) and valence (s = —1) 7r-bands, r/ is a convergence factor, and f s . s '(k,q n) 


describes the overlap between the electron and hole wavefunctions (Eq. B2 We use the full 


7r-band wavefunctions (as opposed to linearized band^) to compute the responses of both 


graphene and graphite. All of our calculations use 0 = 0° (Fig. 2b), since the experimental 
data was found to be insensitive to this angle at the low momenum values studied. We 
performed calculations for both 30° and 60° directions as well, and the differences were 
found to be within numerical uncertainty. The results, shown in Fig. 3 for a convergence 
factor of rj = 0.01 eV, agree closely with those of Ref. 16 . 

For graphite, there are four symmetry-inequivalent sublattices, A and B, leading to four 
basis functions rather than the two of graphene. There is no simple analytic expression for 


these wavefunctions, so to perform the integral in Eq. B4 we simply find the eigenfunctions 
of the graphite Hamiltonian, and use those to compute the corresponding energies and 
overlap functions f s . s '. It is then straightforward to numerically integrate Eq. |B4| over the 
3D Brillouin zone of graphite to obtain n 3D (qr, co). 


Appendix C: Dielectric function 

It is important in our analysis to properly describe the interaction between screening 
processes involving 7r electrons and cr electrons. To do so, we treat the electrons from 
the 7r-bands and cr-bands as two separate fermionic species. In general, we can write the 
Hamiltonian for such a system as 


H = H 0 + Hj + H v , 


(Cl) 
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where 


Ho 'y ^ -SfeCfco-Cfco- + y ] £k9kcr9k<T (^ 2 ) 

ka=^i ka=ti 

describes the noninteracting Hamiltonian for the fermionic operators of the n- and o- 
electrons, g(g^) and c(c7), respectively, with e k and E k their corresponding energy spectra. 

Hi =\ V (9 )M“Q)pM + V (9)M-q)Mq) 


^2v(q)p*(-q)pAq), 


(C3) 


is the Coulomb interaction energy, with p n and p a the corresponding particle density oper¬ 
ators of the 7T and a electrons, V(q) is the Coulomb interaction (see main text), and 


Hu = ^2u{q) [pM + pAq)} 


(C4) 


is an external potential that couples to the density. 

At the RPA level, the equations of motion for the density of 7 r and a electrons 

■ d(pM) 


dt 

, d(pq(q)) 

dt 


= {[H,pA q)]) 

= ([h,pA q)]) 


can be calculated in linear response theorj^S and written in matrix form 


as 


(Pa(q,u)) 

(pn(<l,v)) 


u ( q )_ M I 


where 


and 


M = 


\ UAq,u), 

1 - v(q)nAq, w) V{q)Ua{q, u) 


V{q)Tl n (q,u) 


V(q)UAq,ui) 


(C5) 

(C6) 

(C7) 

(C8) 


e(q, oj) = det M 

= 1 - V{q)UAq,uj) -V(q)U a (q,^) 

is the dielectric function of the system, where I R and 1 b are the RPA polarization functions 


of the 7T- and cr-electrons, respectively. Equation (C7) is completely general and is valid for 


any two band system. The dielectric function can be equivalently written as 


e(q,u) = K a (q) “ V{q)UAq, w), 
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where we have defined 


M<?) = 1 - V{q)n a (q,uj) 


(C9) 


as the background dielectric constant due to the a-electrons. In this way we recover Eq. [7] 
for the total dielectric function of graphene and graphite. 


Eq. C7 simplifies in two equations: 


and 


(p a (q,iuj)) = U(q) 


= U(q) 


n a(q,uj) 

e(q,u) 

n „(q,u) 


(CIO) 


(Cll) 


e{q,u) 

The first equation defines the susceptibility of the sigma electrons in the presence of screening 
effects from both bands, 


X<r(q,u) = 


U a (q,u) 


(C12) 


e(q,iu) ' 

The second one gives the susceptibility of the n electrons in the presence of the sigma bands: 


xA<i,u) 


n n(q,uj) 

e(q,u) 


(C13) 


The full susceptibility is defined by the response function including both cr and 7r electrons 
density fluctuations, 

(MQ, H) + (pAq, iu)) = U(q)x(q, w), (C14) 


where 


x{q, A = x*{q> u ) + Xa(q, A 


n njq.uj) + n g(q,uj) 

e(q,u) 


(C15) 
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FIG. 1. The imaginary part of the experimental density response of graphite and graphene com¬ 
pared to results from RPA calculations, (a) — Im% 3 £>(q,uj) from IXS measurements (blue points), 
RPA calculations including interlayer hopping (green line) and omitting interlayer hopping (red 
line), plotted against energy for selected in-plane momenta, q. The RPA curves had to be multi¬ 
plied by an arbitrary scale factor to allow visual comparison to the experiment. The black, vertical 
lines indicate the energy, hvpq, at which the onset of the particle-hole continuum is expected 
in RPA. (b) — Im% 2 d(q,w) determined from IXS measurements from Eq. 3 (blue points), RPA 
calculations including interlayer hopping (red line) and omitting interlayer hopping (orange line), 
plotted against energy for selected in-plane momenta, q. The black line, again, indicates the on¬ 
set energy of the continuum expected from RPA. Experimental error bars are derived from both 
Poisson statistics and our fitting parameters (see Appendix [A]). 
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FIG. 2 . | (a) | Crystal structure of ABA-stacked graphite defining the hopping parameters t, 71 , and 
73 . |(b) | Brillouin zone of graphene defining the in-plane azimuthal angle of the basal plane. The 
spectra were found to be independent of this angle in the range of momenta measured in this study. 
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FIG. 3. Imaginary part of the polarization function of graphite, — Imll 3 £)(q, u>), from experiment 
(blue points), compared to the results of RPA calculations both with (green line) and without (red 
line) interlayer hopping, 7 ^ 3 . Within the assumptions described in Section III, Imll 3 £>(g,u;) should 
be the same as ImIl 2 _D(q,C(;) apart from an overall factor of the interlayer spacing, d. Note that 
the magnitude discrepancy observed in Fig. 1 is not observed here, because polarization functions 
from different excitations combine additively (Eq. 7) 
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FIG. 4. |(a)| Asymptotic screening and the effective fine structure constant of graphene, (a) Com¬ 
bined plot of the static ratios — lmx3D(q,0)/q and —lm.X2D{q,0)/qd from experiment and from 
RPA calculations both with ( 71,3 7 ^ 0) and without ( 71,3 = 0) interlayer hopping. The results 
from our previous study 1 are also plotted, for comparison. This plot shows that all quantities 
exhibit the proper asymptotic properties in the regime of small q. (b) Momentum dependence of 
the magnitude of the static value of the fine structure constant of graphene. This plot suggests a 
limiting value in the range 0.25-0.35. 
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(c) (d) 

FIG. 5. Illustration of data processing steps described in Appendix A. (a) Subtraction of the 
elastic scattering from the IXS spectra, (red points) Raw IXS data for q = 0.212A -1 . (blue line) 
Fit to the elastic peak using a pseudo-Voigt function, (black points) Resulting spectrum, showing 
the onset of n —> tt* excitations at ~ 1 eV. (b) Matching of this spectrum to wide-range data 
from our previous studj®, to aid in evaluting the sum rule integral, (c) Same illustration as in (a) 
for momentum q = 0.282A -1 . (d) Matching of the q = 0.282A -1 spectrum to data from Ref.^, 
illustrating the mismatch present in a few of the spectra. This mismatch is incorporated into the 
determination of the error bars on the quantity Imy 3 £)(g,u;). 
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